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Abstract 
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Xbj J/ipJ/^ in the framework of nonrelativistic QCD factorization. The relativistic effect 
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is predicted to be of order 10 -5 for X60 62 — > J/fpJ/' l P, but 10~ n for Xbi — > J/ipJ/ip, since there 
9-; is only electromagnetism contribution in this channel. We predict it is possible to observe these 

Q-[ reactions in LHC. Finally, we estimate the decay width and branching ratio of Xcj 

in the constituent quark model by our formula at the leading-order of relativistic correction and 
electromagnetic correction. The obtained branch ratio of Xc0,2 is in agreement with the 

experimental measurement in the order of magnitude. 
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I. INTRODUCTION 



The nonrelativistic nature of heavy quarkonium provides people with a valid component 
to understand the nonrelativistic effect of QCD. Since the quarkonium includes several 
well-separated scales which contain both hard scale and soft scale, the studies of hadronic 
reactions involving heavy quarkonium are significative of giving us an insight into both 
perturbative and nonperturbative QCD. 

The running of B factories with high luminosity has made the measurement of hadronic- 
exclusive processes feasible and given the results that raised new challenges to the existing 
theory. One of the largest puzzles is the cross section for exclusive double charmonium- 
production in annihilation i.e. e + e~ — > J /if) + r] c , at the B factory energy of \fs = 
10.6 GeV(Ref. [lH3j]) is about an order of magnitude larger than the leading-order non- 
relativistic QCD(NRQCD) predictions in Refs. 

The NRQCD factorization [3] as an outstanding effective field theory approach to dealing 
with the physical problem involving multiple scales provides us with a systematic framework 
to deal with the exclusive quarkonium production process; the amplitude can be factorized 
as the products of the short- distance, perturbatively calculable coefficients and the long- 
distance, nonperturbative NRQCD matrix elements which are universal for all processes. 
The short- distance coefficient can be expanded by the order in a s , while the NRQCD matrix 
elements are organized as a series in the relative velocity v of the heavy quark, so one can 
improve the NRQCD predictions simultaneously in a s and v. Considering the next-to- 
leading-order perturbative corrections to e + e~ — > Jhb + rj c , the discrepancy between theory 
and experiment measurement is greatly alleviated [8-131 . 



In recent years, the studies of double-charmonium production at B factory have already 
obtained huge theoretical progress in understanding hadronic exclusive processes with char- 
monium production. These facts inspire us to explore more analogous and valuable processes 
to add to our knowledge about these exclusive processes involving quarkonium. Similar 
to e + e~ annihilation, double-charmonium production in bottomonium decays can also be 
used to study the dynamics of hard exclusive processes and the structure of charmonium 
mesons; the reason is that the bottomonium can decay into not only light hadron but 
also charmonium because its mass is heavier than charmonium and close to the energy 
of B factory. Although the investigations about bottomonium decay are not as intense as 
double-charmonium production in annihilation, there are several works focusing on these de- 
cays (l~3-[22|. In Ref. [13-22], the double-charmonium production in exclusive bottomonium 
decays rjb — > J/ipJ/ip has been investigated. These researches groped the potential of the 
discovery of this hadronic decay channel in experimentation. As the 5-wave bottomonium 
decays into double J /if) have been researched, it is natural to carry on a further study of the 
similar P-wave bottomonium decays into double J /if). 

This paper is intended to deal with the P-wave bottomonium decay process Xbj 
J /if) J /if>, which is a part of the processes dealt with in Refs. 14|, 231. In Ref.[I3], light-cone 
method is employed, and only the QCD process is considered. However, using NRQCD 
factorization to handle this type of process with the initial and final states all involving 
heavy quarkonium is more natural. In this paper, we calculate the relativistic corrections 
to this process in the NRQCD frame; the pure QED process is also considered. In fact, this 
work is an extended and updated version of Ref . j24j . 

This paper is organized as follows. In Sec. HI1 we describe the NRQCD factorization 
formula relevant to this work and compare our matching scheme with the orthodox doctrine. 
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We also present a detailed description on how to determine the short-distance coefficients 
through relative order v 2 in Xbj J/i/jJ/i/j- In Sec. IHH we perform calculations of the 
amplitudes for Xbj — > J/ijjJ/ij) that include the relativistic corrections to the order of v 2 
of J/if), helicity amplitudes and electromagnetism contributions. In Sec. HV] we apply our 
formulas to investigate the phenomenological impact of QCD and QED corrections to the 
decay width and branching ration of Xbj J/ipJ/ip- By analyzing the numerical results, We 
predict the possibility of the observation of Xbj J/ipJ/4> decay in the experimentation. 
Finally, we summarize our results in Sec. [V] 



II. NRQCD FACTORIZATION AND MATCHING STRATEGY 

Heavy quarkonium is a nonrelativistic system involving multiple scales. Because the 
velocity v of the heavy quark is much less than 1, the following scales are well separated: the 
heavy-quark mass tuq, the relative momentum rriQV, and the binding energy ttiqv 2 . NRQCD 
factorization provides us with a valid effective field theory to separate the scale tuq from the 
others. In the NRQCD factorization formula, the contribution from the scale that is less than 
rriQ is kept in the Lagrange operators, while the contribution from the scale that is more than 
uiq is absorbed in the Wilson coefficients. With NRQCD factorization, the decay width can 
be expressed as the sum over qq channels of products of a long-distance-physics-insensitive 
short-distance coefficient and a process independent long-distance non-perturbative matrix 
element. 

For a Q(p)Q(p) pair with total momentum P and relative momentum q we express the 
momenta of c and c in perturbative matching: 

P P 

p = j + q, p=^-q, (i) 

where q and P satisfy Pi - qi — 0, and p 2 = p 2 = rriq. P is the true total momentum of the 
cc pair, P = p + p, with invariant mass of 2E q . In the rest frame of the cc pair, the explicit 
components of the momenta are P^ = (2E q , 0), q^ — (0, q), p^ — (E q , q), and p^ = (E q , — q), 
respectively. 

To be accurate to order v 2 , there are two methods that can be used to match the short- 
distance coefficients. One is the traditional matching method, in which we need to expand 
the energy of the Q or the Q in the QQ rest frame E g around the pole mass m c in power 
series of q 2 , 



E q = m c + ^- + 0^). (2) 



2m c 

The other method which we used in this paper is to expand every occurrence of m c in the 
amplitude in terms of q 2 /i? 2 , while keeping E q intact: 

m c = E q -^- + 0(^). (3) 

In the first method, when summing the polarization states of cc( 3 Si), bb( 3 Pj), there are new 
factors of E q , which are regenerated in the squared amplitude, we have to reexpand these 
occurring E q factors again and realign the corresponding terms from the leading order (LO) 
matrix element squared to the relativistic correction piece. The second method avoids many 
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complications that emerged in the first one and eliminates the task of matching the cross 
section to the amplitude squared. 

There is another question that needs to be considered. In our matching method, m c has 
been eliminated in favor of E q in the physical matrix element squared, then we need decide 
which value of E q should be taken to give the prediction. After comparing the Eq. (T5]) which 
comes from simple nonrelativistic kinematics and the G-K relation [25| : 

+ + 0(A (4) 

where {v 2 )j/^ is a dimensionless ratio of the vacuum matrix elements defined as follow: 



, 2) (JMA)|^(-i D )W(A)x|0) 
V)J I^ ml{J/i>(\)Wve(X)x\0) * 



We can see that theoretical consistency requires that E q can be fixed in an unambiguous 
manner, i.e. E q appearing everywhere in the short- distance coefficients can be replace by 
Mj/ip/2. By this way, the relativistic effects in phase space integrals are automatically 
incorporated. In addition, choosing Mj/^ as the input parameter is better than m c since the 
mass of J/ip is known rather precisely while the charm quark mass is ambiguously defined. 

Since we no longer need to worry about the complication from the phase space integral and 
sum of the polarization states, we can match the short-distance coefficients at the amplitude 
level. Because they are similar for the Xw,6i,fe2 — > J/ipJ/ip, here we take the Xbo as an 
example to demonstrate how to match the short-distance coefficients. At the leading order 
of a s and the order of v 2 , the amplitude Ai[xbo J/^J/'fp) in terms of the vacuum-to- J/ip 
and x&o _ to-vacuum matrix elements can be written as: 

m,n=l , . x 2m 



M Xb0 = ym^ Q ^2MJi^^2Mj^ 2 J2 c°mn(J/H^ a ■ e(Ai) X |0) 

m,n=0 ^ * 

- l -Dj CT . e (A 2 ) x |0)- 7 =(0| X t (-fB- CT )^| Xfco ), (6) 

where (9 mn are the corresponding short- distance coefficients, which are Lorentz scalars formed 
by various kinematic invariants in the reaction. In particular, they also depend explicitly 
on the helicity A of J/ip. For the Lorentz-invariant amplitude in the left-hand side of 
Eq. ([H]), -M-Xboi ft i s mos t natural to assume relativistic normalization for each particle 
state, since the squared amplitude needs to be folded with the relativistic phase space 
integral to obtain the physical decay width. However, in the right-hand side of Eq. (JH]), 
the J/ip and XbJ state appearing in the NRQCD matrix elements conventionally assume 
the nonrelativistic normalization. To compensate this difference, one must insert a factor 
^2Mj^ v ^M^ v /2M^7 in the right side of Eq. ©. 

To determine the coefficients c^ n we follow the moral that these short- distance coefficients 
are insensitive to the long-distance confinement effects, so one can replace the physical J/ip 
state by a free cc pair of quantum number 3 Si, and replace the physical XbJ state by a 
free bb pair of quantum number 3 Pj, by which the NRQCD operator matrix elements can 
be perturbatively calculated. The short- distance coefficients c J mn can then be read off by 
comparing the QCD amplitude for M[bb( 3 Pj) — > cc( 3 Si, Pi, Ai) + cc( 3 Si, P 2 , A 2 )] and the 
corresponding NRQCD factorization formula. 
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After simple NRQCD calculation, the perturbative NRQCD matrix elements for QQ( 3 S\) 
and QQ( 3 Pi) states at leading order are simply expressed as: 



(QQ( 3 S 1 )\^e-(Tx\0} = V^Nc, 

(QQ( 3 Si)\^ (-1%) <>■ ■ e *|o> = v^q 2 

1 

71 



<0| X T (-f£> . (j W\QQ(*P Q )) = ^2N C |q 



(7a) 
(7b) 

(7c) 



where y2N c comes from the spin and color factors of the NRQCD matrix elements and the 
state IQQ^^Lj)) is nonrelativistically normalized. 

Using the formula above, we can obtain the partonic level amplitude A4[bb( 3 Po) — > 
cc( 3 Si, Pi, Ai) + cc( 3 Si, P2, X 2 )} expanded to the order v 2 : 



1 1 , 

M, Po = 8E( qi )E(q 2 )E(Q)J2T, c °rnn^( Q \xH-lD-a)i;\bK 3 P )) 



(8a) 



cc( 3 S l , X 1 )\^ [ --D ) a- ■ e(X 1 ) X \0)(cc( 3 S 1 , A 2 )|-0 T ( ) a ■ e(A 2 ) X |0) 



m=0 n=0 
, 2m 



I tt 



2n 



(2N c f 2 8E( qi )E{q 2 )E(Q) 



, of , q| , 

C 00 + C 10 o+ C 01 m 2 + 



(8b) 



In Eq. (ED, we use relativistic normalization for the cc and 66 states in the computation 
of the QCD amplitude and nonrelativistic normalization in the NRQCD matrix elements. 
Consequently, a factor 2E q appears in the second expression of Eq. (jSJ). 

From Eq. (jSJ), it is straightforward to extract the short-distance coefficients c^ n : 



in 



2{m+n) 



Qm gn 



m\n\ deft" 1 dq 2 " 



M \bb{ 3 Pj) -> ccf^, P 1; AQ + cc( 3 S u P 2 , A 2 )] 
(2N c ) 3 / 2 8E( qi )E(q 2 )E(Q) 



(9) 



q?=q|=o 



We can derive the LO coefficient Cq by putting q — > in the amplitude and equating E q 
and m c . While deducing the coefficient Cq^c/q, we need first expand the amplitude to the 
first order in q 2 prior to taking the q — > limit. 



III. COLOR-SINGLET MODEL CALCULATION 



In this section, we present a calculation for Xbj — > J/^J/^ in perturbative QCD scheme. 
As we discussed in Sec. HI1 the short- distance coefficients are insensitive to the long-distance 
confinement effects, and to obtain the short-distance coefficients Cq and c^, cf , we need 
only to compare the QCD amplitude and the corresponding NRQCD factorization formula 
in QQ state. So we replace the physical J/ip, Xbj states by a free cc pair of quantum 
number 3 Si and a 66 pair of quantum number 3 Pj, respectively, and we compute the QQ 
analog amplitude Ms Pj of the hadronic amplitude M. XbJ) where the amplitude for QQ level 
perturbative process A^[66( 3 Pj) — > cc( 3 Si,Pi) + cc( 3 Si,P2)] has been aliased as A^3 P/ and 
the hadronic level amplitude Ai[xbj —tJ/ift + J/i/)] has been aliased as M. XbJ - 

The Feynman diagrams for the exclusive process Xbj(Q,q, A) — > J/ip(P 1 , qi, Ai) + 
J/ip(P 2 , q 2 , A 2 ) are shown in Fig. 1. 
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FIG. 1: Leading-order QCD Feynman diagrams that contribute to Xbj ~ > J/^J/^ 



A. order-v 2 QCD amplitude 

The QQ analog A4q( 2s+1 Lj) of the hadronic can be obtained by restricting the QQ to 
have an appropriate spectroscopic state. A given spin state of the color-singlet QQ pair can 
be projected out by replacing u(p)v(p) or v(p)u(p) with a projection matrix that can project 
to a particular spin and color channel. In our case, the spins of J/tp and Xbj are all equal 
to 1, so the projection matrix can be expressed as [4J, Il2 ] 



4y/2E q (E q + m Q ) VK 



^v.jt^gg^a^.^ (10b) 

where P 2 = P 2 /4 = wJq — q 2 , N c = 3, and 1 is the unit color matrix, e is a spin polarization 
vector satisfying P ■ e = and e • e* = —1. 

After projecting out the S'-wave color-singlet spin-triplet of cc and the P-wave color- 
singlet spin-triplet of bb , we can expand A^[66( 3 Pj) — > cc( 3 S\, P%, Ai) + cc( 3 Si, P 2 , A 2 )] 
to the order of v 2 and v 2 ,, where v 2 = q 2 /m 2 . Then we can project out the diagonal, 
antisymmetric and symmetric traceless components of bb( 3 Pj) for J = 0, 1, 2 as 0]. It is 
straightforward to obtain the QQ analog M.ipj as follow: 

iq 4 2 6 (N 2 - 1)M 2 

^p = \^ C NI ,2 W m [*Ko< ■ ^ ~ (v? + vl) [(M 2 bo - 2M 2 m )e\ ■ e* 2 - 5P, • e* 2 P 2 -el]], 

Ms Pl = 0, (11) 
ia 4 2 n / 2 (l - N 2 )e 

Ms P2 = 9s a 1-3/2 iirf ^ [ 3M x, 2 ( 2P i ■ + 2P 2 ■ $Pi% ~ 2P i Pp 2 a e*i ■ 4 ~ M 2 b2 ere?) 

• >A ' M Xb2 

+ (v 2 + v 2 )2M 2 ^ (2PfP 2 CT e* ■ e* 2 + 2M$ M <? x 'e? - 3e?PfP 2 ■ e\ - 3P 1 ■ e^P^)] , 

In Sec. HIl we give the amplitude in hadronic level A4[xbj J/i> + J/i)] and the QQ 
analog amplitude [66( 3 Pj) — > cc( 3 Si, Pi, Ai) +cc( 3 Si, P 2 , A 2 )], respectively, in Eqs. and 
OH]). Then, we can derive the following relationship from the these two equations, taking xwi 
for example: 



M 



XbO 



) 2%{Oi)jM / 2M J/ ^(0 1 ) J/fe / 2M XW (0 1 ) XW 
2A c (2P( gi )) 2 y 2A c (2P(g 2 )) 2 y 2N C (2E(Q)) 2 3p ° 1 ] 
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(Oi)j/* = |(J/^|^ + ^|0)| 2 = §^(0) (13) 
Wx* = ^|(0|x + (-^^)|x W )| 2 = ^l(0) (14) 

As we have discussed, for the precision of our work, we have fixed 2E(q 1 ) = Mj/^, 2E(q 2 ) = 
Mj/^ 2 and 2E{Q) = M XbJ . Again the mass of J/ip is identic, so we have Mj/^ = Mj/^ 2 = 

Substituting Eq. lfiT]) into Eq.( 1T2"j) , we reach the order v 2 QCD amplitude M. XbJ : 

M s Xb0 = ^=M Xw M 2 m {-QM 2 XbQ el ■ el + (v 2 + v 2 ) [(M 2 w - 2M 2 m )e\ ■ e\ - 5P, ■ e* 2 P 2 -el}}, 
M*=0, (15) 



iA 2 y/2. 

1 

6 

+ (v 2 + v 2 )2Mj^ (2Pf P°e\ ■ e* 2 + 2M 2 Xb2 el P e? - 3^PfP 2 ■ e{ - 3P, ■ e^/P^ } 



M Xb2 = -^M Xb2 e pa {3^(2^ ■ e* 2 P 2 p er + 2P 2 ■ e\P^ - 2P[P°e\ ■ e* 2 - M^eTe?) 



where 



Ao = 9 f (iVc Jj (JMI^V . e (A 1 )x|0)(^MI^ t ^ • 6(A 2 )x|0)^=(0|x f (-f D ■ a )^\ Xm ) (16) 



'j/V iW Xbo 



A 2 = £ 2 "} N \ Jl {J/^a ■ e(X l ) X \0)(J/H^o- • e(A 2 ) X |0) ^(0| x t(_iBM6«(A))Vix 62 } 



n^Mj^mZI 2 



'J 



B. Helicity Amplitude 

The polarized decay width and branching ratios can offer more useful information for both 
experimentation and theory, which are lost in the unpolarized ones. As we have gotten the 
amplitude of M.[xbj ~^ J/il>J/ilA, for J = 0, 1, 2, we can easily know the helicity amplitude 
by helicity amplitude formalism [26| . According to Ref . [27[ , we can obtain the corresponding 
helicity amplitude: 

Mg^M^e^dW^e), (17) 

where Ai, A 2 is the helicity of J /if), -MaiA 2 is the reduced helicity amplitude which is a 
function of J, Ai, A 2 and particle masses but is independent of 6 and if, /i is the Xbj spin 
projection on fixed axe, and 9 and (p are polar and azimuthal angles of one of the final J/ip 
in the Xbj rest frame. We can obtained the reduced helicity amplitudes Ai as follows: 
Forxbo, 

M^oixm) = ^ M *>o [12M 2 J2M^ - M 2 J + (v 2 + v 2 )(8Mj^ - 3M x 4 w + 12M 2 ^M 2 J] , 
-^titeo) = ^%M XM M 2 ^ [6M 2 w + (v 2 + v 2 )(2M]^ - M 2 J] , (18) 
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For Xb2, 



MUxm) = M x. 2 [3< 2 (4M^ + MlJ - (v? + v 2 )2M 2 ^(3M 2 M + 4M^)] . 
MU( Xb 2) = 0M X62 [6M^M^-(v 2 1 + vi)M^(M^ 2 +4M^)] , 
M'i,o(X*) = ^Mjf+M^ [l2M* b2 - (v 2 + v^)(M x 2 m + 12M^)] , 

MUM = —J— M L WL - <v? + v 2 )M 2 ,J , (19) 

According to parity invariance, there are only two independent helicity amplitudes for 
Xbo and four independent helicity amplitudes for Xb2- The unpolarized amplitude squared 
can be obtained by integrating the polar angle and summing all the helicity states. The 
relationship between the unpolarized amplitude squared and the helicity amplitude squared 
is: 

\ M xJ = \M s ,o(Xbo)\ 2 + 2\MUxbo)\ 2 (20a) 
\ M xJ = l^,o(X62)| 2 + 2|^t 1 (X62)| 2 + 4|^tote2)| 2 + 2|^t_ 1 (x b2 )| 2 (20b) 

Here we give some general prediction of these helicity amplitudes by hadron helicity 
selection rule. In this exclusive decay process, when we consider the lowest-order strong 
interaction in the limit mb — > oo with m c fixed, we can obtain the asymptotic behavior for 

/ 2\ 2+|A 1+ A 2 | 

Br str [ Xw J/^(A0 + J/V(A 2 )] ~ °?A ( -§ , (21) 

where v c is the relative velocity of cc in J/ - ?/', and the factor v 6 comes from the wave function 
at the origin of J/ip- 

From Eq. ( 12T|) . we can see clearly that the helicity configurations of two J/ip decide the 
scaling behavior of the branching ratio. We figure out that when the decay configuration 
is Ai + A 2 = 0, i.e. (Ai, A 2 ) = (0, 0), and (Ai, A 2 ) = (l, -1), the branching ration exhibits the 
slowest asymptotic decrease Br str ~ l/mf- And we can also expect Br str ~ \fm\ when (Ai, 
A 2 ) = (0, 1) and the Br str ~ 1/mf when (Ai, A 2 ) = (l, 1). In fact, due to the nonzero charm 
mass helicity conservation is violated. 

In particular, the helicity state (0,0) of the channel Xbi J/ipJ/ip is zero because it is 
strictly forbidden due to the conflict between parity and ang ular momentum conservation. 
We call this kind of process an unnatural decay process |28j. 



C. Electromagnetic amplitude 

For completeness, we also consider the electromagnetic contributions to Xbj J/4>J/ip- 
There are four QED diagrams, two of which have the same topology as Fig.l, but with gluons 
replaced by photons, lead to the amplitude that has the same form as Eq. ffl6|) except a 2 is 
replaced by e 2 e 2 a 2 . But their contributions are much more suppressed than those from the 
fragmentation diagrams in Fig. 2, and we will not consider them. 
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We can see that in Fig.2, both J/ip decay from the corresponding virtual photon, so we 
can use their decay constant instead of the vacuum matrix element to describe them. With 
the definition given in Ref. [291 ]. 

<0| c^c\J/4>(P, e)> = f m M m & (22) 

where P is the momentum of J/ip, e is its polarization vector, and fj/$ is the so-called decay 
constant. 




FIG. 2: Lowest-order QED diagrams that contributes to Xbj J/tp + J/ip- Only the 
fragmentation-type diagrams are retained, whereas the other two, which can be obtained by re- 
placing the gluons in Fig.l with photons, have been suppressed. 



We can infer the QED fragmentation contribution to the amplitude in Fig. [2j in the 
following results, we have already included all order relativistic corrections by using the 
decay constant: 

'9 R M 2 

MZ = [( 3M L - " 4Ko + 2(2M^ - 3Ml o ) Pl ■ e* 2 P 2 ■ e*] (23a) 



M Z = m 1 ^ L e a p pa Q^(P 1 ^t P P 1 ■ 4 + P[t\°P 2 ■ e x ) , (23b) 

Xbi 

MZ = i2B 2 V2Mj^M Xb2 e pa (23c) 
[2e\^P 1 ■ e* + 2e* 2 p P[P 2 ■ e{ + (2M^ - m 2 ) e\ p e\° + 2P P P?e\ ■ e* 2 ] , 



where 



Bo = -ele^f)^— J J =T^ L ^\ X \-^D • a )V|X60> , (24a) 

V M Xbj M j/-4> » 



Bi = -eg^eVjM^ ^Tl^^ ^(0|x f (-iB x <t ■ e(A))^| XM > , (24b) 



zjYjVf 2 _ 2jVf 2 )~ 2 

5 2 = -^e 2 e%, E^I^C-^^^W^lXM) • (24c) 



o 
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Following are the helicity amplitudes: 

16 ~ 4 

M™{ X m) = iB ^M XM M^,^-( X6 o)=^o^M Xw M^(8M^-3M^ o ) ) (25a) 

Ml m ( Xb i) = iB x AM^{AM^ - MJJ , (25b) 

16 ~ 8 

Ml m ,{ Xb 2) = iB 2 -=M^M XM , MRixn) = iB 2 —M]^M Xb2 , (25c) 

Mr%{xm) = iB 2 4M* b M 3 j^ , M^ixm) = zB 2 AV2M^M Xb2 (M 2 M - 2M 2 /J , 

The relationship between the unpolarized amplitude squared and the helicity amplitude 
squared is similar to the relation in the QCD section except a new equation for xti- 

\MZ\ 2 = mt m Q ( Xb i)\ 2 (26) 

D. Decay width 

As we have known the amplitude squared, we can easily derive the decay width of Xbj 
J/ipJ/ip. The polarized decay width can be expressed as 

1 1 ^ f 1 dcosO. , .(j) l9 
a.a, = >< o >< $ 2 / -^\M^ iX2 J (27) 



2M X6J 2 V-i 2 
The unpolarized decay width can be expressed as 



^ 1 1 ^ Z" 1 dcosQ s . , .fji l9 

r unp = x - x $ 2 / > LMr\ , 2 28 

p 2M y „ 2 2 / l 2 ^ 1 AlA2 ^' v 7 

Xw J 1 AiA 2 



where = \M^ S ~ M^ em \\ ^ = hi 1 ~ ^ 



XbJ 



IV. NUMERICAL ANALYSIS 
A. Input parameters 

In this section, we will apply the results obtained in Sec. IHII to give some phenomenological 
predictions. Before we carry out the numerical calculations, we need to fix several input 
parameters, such as M XbJ , M J/i} , a, a s , {Ox)j/y, (Oi) XbJ , fj/^p and (v 2 ) J/i} . 

First, we need to fix the values of coupling constants; in our work, we set the running QCD 
strongcoupling constant a s (M XbJ /2) = 0.22 by using the two-loop formula with A^g = 0.338 
GeV ]!, [9| and the electromagnetic fine structure constant a = 1/137. 

Next, The NRQCD matrix element (0\) XbJ can be obtained from the derivative of radial 
wave functions near the origin in the potential model, for P-wave is (Oi) XbJ ~ ^^|-R'ip(0)| 2 
as in Ref. [i, 30] . In Ref. j3l[ the values of \R[ P (0) | 2 for four potentials have been listed; we use 



the value of the Cornell potential, and then obtain {Oi) XbJ = 2.96 GeV . The color-singlet 
NRQCD matrix elements of J/if>, i.e.(Oi)j/$, can be obtained from the electromagnetic 
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decay rate of the J/ip in which the a s leading order is considered. Using the measured 
dielectric width 5.5 KeV, we obtained (Ci)j/^, = 0.268 GeV 3 . 

The values for the physical masses of the involving hadrons are taken from Ref. [32| as 

9.912 GeV^M J/v , = 3.097 GeV. The decay 



M Xb0 = 9.859 GeV, M Xbl = 9.893 GeV, M Xb2 



0.406 GeV as in Ref. 29 



constant for J /ip is taken as fj/^ 

Considering the order v 2 relativistic corrections are calculated, we also need to know the 
value of (v 2 )j/^. Here we adopt the value (v 2 )j/^ = jy 2 )j/^ 2 = 0.225 extracted from the 



recent Cornell-potential-model-based analysis in Ref. [33 



B. The decay width of Xbj J/ipJ/ip 



With the input parameters fixed, we present the numerical results of the polarized and 
unpolarized decay widths for Xbj — > J/tpJ/ip in Table HI 



TABLE I: The polarized and unpolarized decay widths for Xbj J/i^J/tp, where we list the 
leading-order QCD contributions of T(xbj —> J/ipJ/ip), QED corrections, i> 2 -order QCD corrections 
(including LO contributions) and both the QED and QCD corrections. The values are all in units 



of eV. 







To, o 




Ti, o 


ri-i 


^unp 




X60 


4.334 


0.262 






4.859 


LO 


Xbi 














Xb2 


1.251 


0.099 


0.758 


3.881 


12.240 




Xbo 


4.338 


0.266 






4.870 


QED 


Xbi 






2.26 x 10~ 7 




9.04 x 10~ 7 




Xb2 


1.250 


0.098 


0.757 


3.872 


12.216 




Xbo 


5.067 


0.231 






5.530 


v 2 -QCD 


Xbi 














Xb2 


1.078 


0.079 


0.639 


3.439 


10.673 




Xbo 


5.070 


0.235 






5.541 


QED&QCD 


Xbi 






2.26 x 10~ 7 




9.04 x 10~ 7 




Xb2 


1.077 


0.079 


0.638 


3.431 


10.651 



We can see that the relative order-f 2 contributions are very prominent. It is found to 
increase the lower-order prediction for the decay width by about 13.8% for xw — > J/ipJ/ip, 
and about 12.8% for Xb2 J/ipJ/ip- The QED contributions are very small, only about 
0.2% for Xbo -> J/ip J/ip and Xb2 -> J/ipJ/ip- 

We can also see in Table |T] that the polarized decay width is consistent with the helicity 
selection rule in general. As we analyzed in Sec. IHIt the configurations of two J/ip with (Ax, 
^2) = (0, 0) and (Ai, A 2 ) = (l, -1) bear the smallest suppression and contribute the largest 
component to total decay width. The (Ai, A 2 ) = (0, 1) channel gives the second largest 
contribution and the contribution from the (Ai, A2) = (l, 1) channel is the smallest one. So 
we speculate the most observable channels are (Ai, A2) = (l, -1) for Xb2 and (Ai, A2) = (0, 0) 
for Xbo respectively. 
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C. The decay branching ratio of Xbj — > J/iftJ/iJj 



Nevertheless, what we are more interested in is the branching ratio, so we need to know 
the total decay width. According to Ref. [tJ, the decay rate of Xbj into light hadrons (LH) 
is 



r( XbJ -> LH) 



2Im/ 1 ( 3 P 7 



ml 



■{O x {*Pj)) XhJ + 



2Im/ 8 ( 3 50 



mi 



(OsCSo)) 



Xbj 



(29) 



The imaginary parts of /i( 3 Pj) and fs(}So) have been calculated up to order a 3 in 
Ref. 341 ] . but for /i( 3 P ) an d /i( 3 P2) different results were given in Ref.|35j], so we take the 
leading-order results as 



Ftot(.Xbo) 
Ftot(Xbl) 

r*ot(Xfe2) 



64^(0^ Pj)) XbJ A n f Ml (O 8 CS )) XbJ 



3Mt 

XbO 



647ra 2 (0 1 ( 3 P J )) 



Xbj 



3M4 

Xbl 



+ 



16(O^Pj)) XbJ 
n f Ml bi (O^S )) XbJ 
l6{O^Pj)) XbJ 

64 7 ra 2 (0 1 ( 3 P J )> x , 7 (4_ , n f MlJO^S,)) XbJ 

15 



\Q(O^Pj)) XbJ 



(30a) 
(30b) 
(30c) 



0.0021 GeV~ 2 , then we can 



where n f = 4. According to Ref. H, we take ffp"))*" 
obtain the value of total decay. 

Another method is to calculate the decay width of El transitions of Xbj 7T, according 
to the branching ratio of Xbj — > 7^ from Ref. 36[, we can obtain the total decay width of 
Xbj- The El transitions width is defined as j37]~ 



r(l 3 p J ^>l 3 ,S 1 + 7 ) = ^-A; 3 



r 2 (irPn(r)Pio(r) 



(31) 



where is the photon momentum, and P n ;( r ) * s radial wave-function. 

The total decay widths of Xbj obtained from the two methods are listed in Table |TTJ 



TABLE II: Total decay width of Xbj- One is obtained by El transitions, which is from Ref. 37] 
and PDG data. Another is from the Xbj decay into LH in leading order. The values are all in units 
of KeV. 





El 


LH 


r(xfco) 


> 319 


1068 


r(x6i) 


69 


52 


r(x fe2 ) 


124 


317 



From all these values we can obtain the branching ratio in Table IIHI For comparison, we 
juxtapose the decay width of Xb.i —> J/ifrJ/ip from our calculation and the results obtained 
in [23J within NRQCD, where the first errors come from the uncertainties in the distribution 
for final-state mesons, the second in a s , and the last errors are associated with power-law 
corrections. 
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TABLE III: Total decay width and branching ratio of Xbj J/ipJ/ip. Here we give the total 
decay width including QED and order-f 2 relativistic corrections. For comparison, we list the 



results obtained in 
into light hadrons and El transition. 



23J within NRQCD. We also show the branching ratios obtained by Xbj decay 



M 




T(eV) 


r[23](eV) 


Br E1 (10- 5 ) 


Br LH (10- 5 ) 


Xbo ~ 
Xbi - 

Xb2 - 


■»• J/ipJ/^ 


5.54 
9.04 x 10~ 7 
10.6 


27±| 5 ± 19 ± 13 
65±^ ± 46 ± 32 


< 1.7 
1.4 x 10~ 6 
8.6 


0.5 
1.9 x 10~ 6 
3.4 



We can see that the results obtained in 



23| are in general larger than ours. We have 



compared our analytical expressions with 23| and found a difference by a factor of 2 in the 



formula of decay rate which will reduce the result of 23[ half. In addition, the various input 
parameters such as NRQCD matrix elements, a s , and particle masses can bring the large 
uncertainties. Finally, because we expand the pole mass m c around the physical mass Mj/^ 
instead of the traditional matching method, this will decrease our predictions comparing to 
the orthodox method. Although the results seem to have a sensible difference, we believe 
the order of magnitude is correct. 



D. Some predictions of observation potential for Xbj ~^ J/^J/tp 

We use the branching ratio given in the Table II III to explore the possibility for these 
processes to be observed in experiments. Considering the branching ratio for each of the 
decays J/ip — > fi + fi~ is about 6%, we find the branching ratio of Xbj decay to muon pair is 
Br[ X bo -»■ J/ipJ/1/) M ~ 0.6 x 10~ 7 and Br[x b2 -> J/i>J/i/> -» M « 3.1 x 1(T 7 . The 
total cross section for Xbo and Xb2 production at Tevatron energy according to [14J is 

cr(pp Xbo + X) = 250 nb, 

<r(pp -> Xb2 + X) = 320 nb, (32) 

We estimate there are about 50 ~ 150 for Xbo an d 400 ~ 1000 for Xb2 produced events in 
Tevatron Run 2 that achieve an integrated luminosity of about 10 fb -1 by April 2011. 

Similarly, we can combine the cross sections at LHC, which are about 6 times larger than 
the corresponding cross sections at Tevatron [14( , to predict the number of produced events 
that may reach between 1500 and 4500 for Xbo arid between 12000 and 30000 for Xb2 with 
the accumulated luminosity 50 fb -1 of LHC in 2010, and with the acceptance and efficiency 
of detector considered, we expect 15 ~ 45 and 120 ~ 300 observed events per year for Xbo 
and Xb2 respectively . 



E. The decay width and branching ratio of Xcj — > VV(V — > ui, <p) 

The constituent quark model, which treats the light mesons as non-relativistic bound 
states, is also frequently invoked as an alternative method for a quick order-of-magnitude 
estimate. In this sense, the preceding formulas derived for Xbj —> J/^J/^ cari be applied 
to describe the decay processes Xcj VV, once we understand that we are working with 
the constituent quark model. 
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TABLE IV: Parameters used for numerical calculation in Xcj <P4> an d Xcj ™ 



My(GeV) a s \R sV (0)| 2 (GeV 3 ) |^J 2 (GeV 5 ) M Xe0 ( GeV) M Xel (GeV) M Xe2 (GeV) 

XcJ^ojoj 0.78265 0.3 0.11 0.075 3.41475 3.51066 3.55620 

Xcj->4><P 1.019455 0.3 0.19 0.075 3.41475 3.51066 3.55620 



TABLE V: Decay widths and Branching ratios of Xcj —> 00 an d Xcj uuj. The second and third 
rows are the results of our calculation. The next two rows are the results obtained in 1J| within 
"^3" model. The next row is the branching fractions obtained by BESIII 39J, where we take the 
results of combined final state. The last row is the PDG [32J results. In Xcj ww, there are only 
our results, BESIII results and PDG results. 



XcO 



Xcl 



Xc2 



r(keV) 
Br(10" 4 ) 
r(keV) [14] 
Br(10~ 4 ) [14] 
Br(10~ 4 ) [22| (BESIII) 
Br(10" 4 ) [32](PDG) 



3.3 
3.2 
2.10 
3.01 
.0± 0.3 ±0.8 
9.2 ± 1.9 



1.9 x 10~ 7 
2.2 x 10~ 6 



4.4 ±0.2 ±0.5 



5.9 
30 
3.38 
21.3 
10.7 ± 0.3 ± 1.2 
14.8 ±2.8 



XcO 



UJUJ 



Xcl ->■ UUJ 



2.2 x 10~ 7 
2.6 x 10~ 6 
6.0 ±0.2 ±0.7 



Xc2 



UJUJ 



r(keV) 
Br(10" 4 ) 

Br(10~ 4 ) [221 (BESIII) 
Br(10" 4 ) [32](PDG) 



2.3 
2.2 

9.5 ± 0.3 ± 1.1 

22 ±7 



3.2 
16 

8.9 ± 0.3 ± 1.1 
19.0 ±6.0 



We have listed the numerical values for parameters in Table. IIVI 

We take Xcj 00 as an representative. By regarding as a strangeonium, we can directly 
use Eq. (j27p . only with some trivial changes of input parameters. We take the constituent 
quark mass m s ~ = 0.5 GeV. The radial wave function at the origin of 0, |-R^(0)| 2 , 

can be extracted analogously from its measured dielectron width of 0.19 GeV 3 . Taking 

M. 



Xcl) 



3.41475 GeV, and the strong coupling constant a s = 0.3, and only considering the 
contribution of the leading-order QCD and QED contributions, we can obtain the decay 
width and branching ration of Xcj ~* 00 in Table |V] 

It is also interesting to consider the similar decay process Xcj wu, Parallel to the 
preceding procedure, we also give the results in Table [V] 

From the Table [V] we can see our predictions of Br[x c0 — > 00] and Br[x c2 — > 00] are 
3.2 x 10 -4 and 30 x 10 -4 resp ectively. Our results are compatible with the results obtained 
in [13] within "03" model 38 not only in decay width but also in branching ratio. Moreover, 
the branching ratio of Xco,2 — > 00 is also compatible with BESIII [39j measurement and 
PDG 32j in the order of magnitude. For the process Xcj ww, We can learn from Table IVl 
that our result of Br[x c0 — > wu] is close to the observation in BESIII [39], while for Br[x c2 — > 
uu}, our result 
BESIII 



39 



is close to the PDG 
and the result in PDG [3_2 



32 



It is noting that both of the observation in 
the branching ratio of Xco ~^ uuj is large than the 



one of Xc2 —> utu, while our prediction is opposite. This discrepancy may be caused by the 
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model approximation. 



V. CONCLUSION 



The major task of this work is to calculate the first-order relativistic corrections and 
electromagnetic corrections of Xbj J/^J/^P m the context of NRQCD factorization. We 
first introduced NRQCD factorization formula, which is particularly suitable for calculating 
the relativistic corrections to quarkonium production and decay processes in the color-singlet 
channel. Then we compared the orthodox matching strategy with our matching method and 
emphasized that our approach avoids the complications that come from squared amplitude 
and phase-space integral. These two methods are equivalent thanks to the Gremm-Kapustin 
relation and we give our matching schema at amplitude level. 

As a phenomenological application of our results, we calculate the decay width and 
branching ratio of Xbj J/ipJ/ip. The conclusion is that the relativistic contributions 
of Xbj J/iiJ/ip are modest, about 13.8% for Xbo — >• J/ipJ/ip, and about 12.8% for 
Xb2 — > J/tpJ/ip. The QED contributions are very small, only about 0.2% for xw J/ipJ/ip 
and Xb2 — > J/ipJ/ip- The more interesting branching ratios are predicted to be of order 
1CT 5 for Xbo J/ipJ/^i Xb2 J/ipJ/ip, but 1CT 11 for Xbi J/ipJ/ip, since there is only 
electromagnetism contribution in this channel. With the result, we predict it is possible to 
observe these processes in LHC. 

Finally, as an exploratory work, we also have a quick order-of-magnitude estimate of the 
decay width and branching ratio of Xcj ww, 00 by our formula in leading-order QCD and 
QED contribution working with the constituent quark model. The predict of Br[x c o,2 — > 00] 
is consistent with the experimental result in the order of magnitude. 
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